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Abstract: By numerically solving the appropriate Boltzmann equations, we study the pro- 
duction of sterile neutrinos in models with low reheating temperatures. We take into account 
the production in oscillations as well as in direct decays and compute the sterile neutrino 
primordial spectrum, the effective number of neutrino species, and the sterile neutrino con- 
tribution to the mass density of the Universe as a function of the mixing and the reheating 
parameters. It is shown that sterile neutrinos with non-negligible mixing angles do not nec- 
essarily lead to N v ~ 4 and that sterile neutrinos may have the right relic density to explain 
the dark matter of the Universe. If dark matter consists of sterile neutrinos produced in oscil- 
lations, X-rays measurements set a strong limit on the reheating temperature, Tr > 7 MeV. 
We also point out that the direct decay opens up a new production mechanism for sterile 
neutrino dark matter where cosmological constraints can be satisfied. 
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1. Introduction 

Sterile neutrinos likely exist. They can easily be incorporated into the standard model and 
provide the simplest explanation for the existence of neutrino masses. The most important 
parameter associated with sterile neutrinos is probably their mass scale. In seesaw models 
[Q], where sterile neutrinos are simply added to the standard model matter fields in order to 
generate light neutrino masses, sterile neutrino masses are free parameters of the Lagrangian, 
whose values are to be experimentally determined. To account for the neutrino masses in- 
ferred from the solar and atmospheric neutrino experiments, at least two sterile neutrinos 
are required but only mild constraints on the masses or mixing of the sterile neutrinos can 
be derived. And theoretical considerations are of no help either, for heavy as well as light 
sterile neutrinos can be motivated on different grounds || ||]. It seems reasonable, then, to 
consider the sterile neutrino mass scale simply as another free parameter subject to present 
experimental constraints. In this paper, we study sterile neutrinos with masses in the eV-keV 
range. 

Sterile neutrinos with keV masses have indeed been proposed as dark matter candidates 
[Q, ||, |6|. In the early Universe, such sterile neutrinos are produced in active-sterile neutrino 
oscillations and never reach thermal equilibrium. Due to their primordial velocity distribu- 
tion, sterile neutrinos damp inhomogeneities on small scales and therefore behave as warm 
dark matter particles. The mass of dark matter sterile neutrinos is constrained from below 
by the observed clustering on small scales of the Lyman-a forest [J?J. Present bounds give 
m s > 10-14 keV [|. Because of its mixing with active neutrinos, the v s may radiatively 
decay (through v s — » v + 7) producing a monoenergetic photon with E~ ~ m s /2. X-rays 
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measurements may therefore be used to constraint or infer the mass of the sterile neutrino. 
Recent bounds, based on observations of the Virgo and Coma clusters and the X-ray back- 
ground, yield m s < 6-10 keV [10, 11, |l!| and are thus in conflict with the Lyman-a forest 
constraint. That means that the minimal mechanism for sterile neutrino dark matter, based 
on active-sterile oscillations, is already ruled out |nj ||, ^j. 

A possible clue regarding the mass scale of the sterile neutrinos is the result of the LSND 
experiment ||13|| . It found evidence of — ► v e conversion, which is being tested by the 



Fermilab MiniBoone experiment [14]. The LSND signal can be explained by the existence of 
light (m s ~ 1 — 10 eV) sterile neutrinos mixed with v e and Vu @ . In the standard cosmological 
model, such sterile neutrinos generate two important problems: i) They give a contribution to 
£l v larger than that suggested by global fits of CMD and LSS data |15|| . ii) They thermalize in 
the early Universe so that N v ~ 4, in possible conflict with big-bang nucleosynthesis bounds 
[|il| . Recently, the MiniBoone experiment presented its first results (lTj which disfavore even 
more the so-called (3+1) schemes JT3[. It seems, nonetheless, that (3+2) schemes are still 



viable [If]. 

The standard cosmological model, however, has not been tested beyond big bang nu- 
cleosynthesis, for T > 1 MeV. Cosmological models with low reheating temperatures, for 
example, offer a natural and viable alternative to the standard paradigm. In fact, various 
scenarios of physics beyond the standard model, including supersymmetry and superstring 
theories, predict the existence of massive particles with long lifetimes that decay about the 
big bang nucleosynthesis epoch, inducing a low reheating temperature and modifying the 
initial conditions of the standard cosmology. Over the years, different issues related to these 
models have been studied in the literature [19, 20, 21]. In this paper we consider the possible 
interplay between sterile neutrinos and models with low reheating temperatures. On the one 
hand, sterile neutrinos may serve as probes of the early Universe and constrain the reheating 
temperature. On the other hand, models with low reheating temperatures may alleviate some 
of the problems associated with sterile neutrinos, suppressing their abundance or modifying 
the standard relation between the sterile neutrino relic density and the mixing parameters. 

So far, a detailed analysis of these effects have not been presented in the literature. 
Cosmologies with low reheating temperatures were suggested, in jl6||, as a possible way to 



accommodate the LSND signal and big bang nucleosynthesis, whereas in [22], several simpli- 
fying assumptions -not all of them justified- were used to obtain and analytic estimation of 
the sterile neutrinos produced in oscillations. In this paper, we numerically solve the equa- 
tions that determine the sterile neutrino distribution function in models with low reheating 
temperatures. Two different sources of sterile neutrinos are taken into account: active-sterile 
oscillations and the direct decay of the field responsible for the reheating process. We compute 
different observables related to the sterile neutrino, including its spectrum and relic density, as 
a function of the reheating parameters and the mixing angle and mass of the sterile neutrino. 

In the next section we describe the reheating process and introduce the different equations 
that are relevant for the production of sterile neutrinos. Then, the behavior of active neutrinos 
in models with low reheating temperatures will be briefly reviewed. In section |3], we study 
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in detail the production of sterile neutrinos as a result of active-sterile neutrino oscillations 
for different mixing and reheating parameters. We show that N v ~ 3 can be obtained even 
for sterile neutrinos with relatively large mixing angles and that dark matter sterile neutrinos 
provide a strong constraint on the reheating temperature. Finally, in section ||, we include 
the production of sterile neutrinos through the direct decay of the scalar field and study 
the resulting sterile neutrino spectrum and relic density. We observe that sterile neutrinos 
produced in decays may account for the dark matter and avoid the Lyman-a and X-ray 
constraints. 



2. The reheating process 

Reheating is defined as the transition period between a Universe dominated by a unstable non- 
relativistic particle, cp, and the radiation dominated Universe. In the standard cosmological 
model reheating is assumed to occur only after inflation, but in general, additional reheating 
phases not related to inflation are also possible and our discussion applies equally to them. 
During reheating the dynamics of the Universe is rather involved. The energy density per 
comoving volume of the non-relativistic particle decreases as e~ r ^* -with the <fi decay 
width- whereas the light decay products of the 4> field thermalize. Their temperature quickly 
reaches a maximum value T max and then decreases as T oc a" 3 / 8 §,asa result of the 
continuous entropy release. During this time the relation between the expansion rate and 
the temperature is neither that of a matter-dominated universe (H oc T 3 / 2 ) nor that of a 
radiation-dominated Universe (H oc T 4 ) but it is given instead by H oc T 4 . Thus, at a given 
temperature the Universe expands faster during reheating than in the radiation-dominated 
era. This unusual behavior continues until t ~ r^ 1 , when the radiation dominated phase 
commences with temperature Tr. From then on, that is for T < Tr, the evolution of the 
Universe proceeds as in the standard scenario but with initial conditions determined by the 
reheating process. 

The success of standard big bang nucleosynthesis provides the strongest constraint on 
Tr. Electrons and photons interact electromagnetically and consequently have large cre- 
ation, annihilation and scattering rates that keep them in equilibrium even during reheating. 
Neutrinos, on the contrary, can interact only through the weak interactions and are slowly 
produced in electron-positron annihilations. Since big bang nucleosynthesis requires a ther- 
mal neutrino spectrum, Tr should be high enough to allow the thermalization of the neutrino 
sea. Given that, in the standard cosmology, neutrinos decouple from the thermal plasma at 
T ~ 2 — 3 MeV, it can be estimated that they will not thermalize if Tr < few MeV. Indeed, 
detailed calculations give T > 2 — 4 MeV 21] as the present bound. In this paper, we 



consider models with reheating temperatures below 10 MeV. 

Let us know formulate the equations that describe the reheating process, and in partic- 
ular, the production of sterile neutrinos at low reheating temperatures. We denote by cp the 
unstable non-relativistic particle that initially dominates the energy density of the Universe. 
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Its energy density, p$, evolves according to 

< ^L = -T m -Wp <j) (2.1) 

where H is the Hubble parameter and is the 4> decay width. 

The energy-momentum conservation equation in the expanding universe is 

^- = -3H(p T + P T ) (2.2) 

where px and Pt denote respectively the total energy density and the total pressure. At the 
low temperatures we allow for, only the scalar field, electrons, photons, and neutrinos are 
present in the plasma. Denoting by p u the energy density in active and sterile neutrinos, we 
have that 

Pr{t) = Pcj> + p 7 + p e + Pv (2.3) 

and an analogous expression holds for Pt- Equation (|2.2| ) can be rewritten as an evolution 
equation for the (photon) temperature as 



cffiy __ + AH Pl + 3H(p e + P e ) + 4Hp„ + dp v /dt 

dt ~ dpj/dT^ + dp e /dT 1 

H, the hubble parameter, is given by the Friedmann equation, 



(2.4) 



with a the scale factor and Mp the Planck mass . 

We follow the evolution of active neutrinos by solving the momentum-dependent Boltz- 
mann equation 

for v e and (f Ur = f Vfl )- C colli the total collision term, describes neutrino annihilations and 
scatterings. The following processes are taken into account in our calculations: 

v% + v% <-> e + + e" (2.7) 
Vi + e ± <r-> Vi + e ± . (2.8) 

The collision terms associated with these processes are complicated, involving nine-dimensional 
integrations over momentum space. But they can be simplified to one-dimensional integrals 
by neglecting m e and assuming that electrons obey the Boltzmann distribution [19|. Since 



the error due to the above approximations is small (less than few percent), we will use the 
one-dimensional form of the collision terms. 

Regarding the sterile neutrinos, we will consider the simplifying limit of two neutrino 
(active-sterile) mixing. That is, we assume one sterile neutrino, v s , that mixes predominantly 
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with a single active flavor v a (a = e,/x, t). In consequence, the transformation between the 
flavor and the mass bases can be written as 

\u a ) = cos9 \v\) + sin# ^2) (2-9) 
\v s ) = -svD.6\vx) +cos0|i/ 2 ) (2.10) 

where \u\) and \v<2) are neutrino mass eigenstates with masses m\ and m2, respectively. 9, the 
mixing angle, parameterizes the magnitude of the mixing between the active and the sterile 
neutrino. For the small mixing angles we deal with, |z/ 2 } practically coincides with |i/ s ), so we 
will use m s instead of m 2 to denote the mass of the eigenstate that is predominantly sterile. 

The sterile neutrino distribution function also follows a Boltzmann equation like (|2[ 
The collision term for v n <^> oscillations is Hi: 



1 T a (p)A 2 (p)sm 2 28 
4 A 2 {p) sin 2 29 + D 2 (p) + [A(p) cos 29 - V T {p)f 



where A(p) = m 2 /2p, T a is the f a total interaction rate, D(p) = T a /2 is the quantum 
damping rate, and V T is the thermal potential. 

In addition to oscillations, we also consider the production of sterile neutrinos through 
the direct decay <f> ~> v s v s . Since <f> is nonrelativistic, each sterile neutrino is born with 
momentum m^/2 and the collision integral becomes 

2vr 2 

where b is the branching ratio into sterile neutrinos, and , are respectively the <f> mass 
and number density. 

As initial conditions we assume that at early times the energy-density of the Universe is 
dominated by (j), and that active and sterile neutrinos are absent from the primordial plasma. 
As long as the maximum temperature reached by the plasma (T max p0| ) is large enough, 
the final outcome is independent of the initial conditions. We found that T max ~ 20 MeV is 
enough to guarantee such independence. 

Our analysis can naturally be divided into two parts: production in oscillations only 
(6 = 0), and production in oscillations and decay (67^ 0). In the first case, to be investigated 
in section ^j, the parameters that enter into the above equations are m s , sin 2 29, and Y^. 
It is customary to trade Y^ with the cosmological parameter Tr -known as the reheating 
temperature- through the relations 

I> = 3H(T R ) (2.13) 

and 

7% (Sn 3 3 ,^ 1/2 



"™ = 3 ^Hfl (2 ' 14) 
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with = 10.75. These equations establish a one-to-one correspondence between and Tr. 
In the second case, when sterile neutrinos are also produced in decays (6^0), the results will 
depend additionally on b and m^. Section || deals with this interesting possibility. 

For a given set of mixing and reheating parameters, we simultaneously follow the evo- 
lution of p^, Try, fu e {jp)i fu^ip), an d fu s (p) from the matter dominated era well into the 
radiation-dominated Universe, until the distribution functions reach their asymptotic values 
(T < 0.1 MeV). The main output from this system of equations are the neutrino distribution 
functions, which can be used to compute several observables. Big bang nucleosynthesis, for 
instance, is sensitive to the relativistic energy density in neutrinos. This quantity is usually 
parameterized in units of the energy density of a standard model neutrino, p Uo , and denoted 
by N u , 

Since sterile neutrinos are dark matter candidates, it is also important to compute their relic 
abundance, 

n s = ^£ (2.16) 

Pc 

where m s , n s are respectively the mass and number density of the sterile neutrinos, and p c is 
the critical density of the Universe. 

3. Active neutrinos and low Tr 

The evolution of the sterile neutrino distribution function strongly depends on the corre- 
sponding function of the active neutrino flavor with which it mixes and it is in many ways 
analogous to it. Before considering sterile neutrinos, it is therefore appropriate to briefly 
review the salient features related to the behavior of active neutrinos in models with low 
reheating temperatures. 

Figure [I] shows the evolution of the electron neutrino number density (normalized to the 
equilibrium density) as a function of the temperature for different reheating temperatures. 
The pattern is clear. At high temperatures, T ^> Tr, neutrinos are out of equilibrium and 
n Ue /n eq continually decreases with time until T ~ Tr is reached. For T < Tr, neutrinos 
evolve as in the radiation dominated but with a non-equilibrium initial condition (n Ue (Tji) ^ 
n e q(Tji)). If Tr is large enough, neutrinos will be able to recover the equilibrium distribution 
before decoupling from the thermal plasma. Such event, illustrated by the line Tr = 8 MeV 
in figure |l], would be indistinguishable from the standard cosmology. For smaller reheating 
temperatures, on the other hand, neutrinos never reach the equilibrium distribution and 
decouple from the plasma with a smaller abundance than in the standard scenario. That is 
exactly what happens, for instance, if Tr < 4 MeV (see figure [l]). Note nonetheless that even 
for Tr = 3 MeV the asymptotic deviation from the standard prediction amounts to less than 
10%. 

Because muons are not present in the thermal plasma at low temperatures, muon neutri- 
nos can only be produced in neutral-current interactions. Consequently, the muon neutrino 
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Figure 1: The evolution of the electron neu- 
trino number density as a function of the pho- 
ton temperature for different reheating tem- 
peratures. 
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Figure 2: The evolution of the muon (or 
tau) neutrino number density as a function of 
the photon temperature for different reheating 
temperatures. 



deviates from equilibrium farther than the electron neutrino, as revealed in figure §. Indeed, 
for Tr = 3 MeV the deviation from the standard prediction amounts to 50%. 

The effects of the reheating process can also be seen in the primordial neutrino spectrum. 
A equilibrium spectrum with T u = T 7 /1.4 is expected in the standard cosmological model. 
Figure |3] shows the Vn primordial energy spectrum for different values of Tr as a function of 
p/T-f. The deviation from equilibrium is clearly visible for the smaller reheating temperatures. 




Figure 3: The primordial energy spectrum of the muon neutrino as a function of pjT 1 for different 
reheating temperatures. 
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Figure 4: The evolution of the sterile neu- 
trino number density as a function of the pho- 
ton temperature for different reheating tem- 
peratures and sin 2 29 = 10~ 2 . 
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Figure 5: The evolution of the sterile neu- 
trino number density as a function of the pho- 
ton temperature for different mixing angles 
and T R = 4 MeV. 



4. Sterile neutrino production in oscillations 

Let us now consider the production of sterile neutrinos through active-sterile neutrino oscil- 
lations. For simplicity we will consider mixing with the electron neutrino only so that sin 2 29 
denotes the mixing angle between v e and v s . We are then left with 3 parameters that deter- 
mine all the observables: Tr, sin 2 29, and m s . In this section we study how these parameters 
affect f Us , N v , and Cl Vs . 

The evolution of the sterile neutrino number density follows a pattern similar to that of 
the active neutrinos. Figure || shows n Us /n eq as a function of the temperature for different 
values of Tr and sin 2 29 = 10 -2 . Sterile neutrinos are always out of equilibrium and n Vs /n eq 
decreases with time during the reheating phase, reaching its minimum value at T ~ Tr. 
At T < Tr, the universe is radiation dominated and the sterile neutrino population slightly 
increases, in part as a result of the corresponding increase in n Ue (see figure [l]). The asymptotic 
value of n u Jn eq , however, differs very little from its value at Tr. 



Note that this result is at odds with 22], where it was assumed that the production of 
sterile neutrinos starts at Tr. Actually, as we have seen, sterile neutrinos are slowly created 
during the (f) dominated era and only a small fraction of them are produced after Tr. 

For the range of sterile neutrino masses considered, n v Jn eq does not depend on m s . 
Thus, the other relevant dependence to investigate is that with sin 2 29. In figure ||, n Ua /n eq 
is shown as a function of the temperature for Tr = 4 MeV and different mixing angles. As 
expected, the smaller the mixing angle the smaller n Us /n eq . Indeed, for small mixing angles 
(sin 2 29 < 10~ 2 ), n Ua /n eq oc sin 2 29, as seen in figure [|. Such proportionality is expected when 
f Us can be neglected with respect to f Ue in equation ( p.ll| ). At large mixing angles f Us may 
become comparable with f Ve and the above relation no longer holds. Neglecting f Vs in ( |2.11 ), 
therefore, is not a good approximation for sterile neutrinos with large mixing angles. 
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The primordial energy spectrum of the sterile neutrino is shown in figures ^ and [7| for 
different values of Tr and sin 2 29. It is certainly non-thermal and is strongly suppressed for 
low reheating temperatures or small mixing angles. 

Standard big bang nucleosynthesis is a powerful cosmological probe of active and sterile 
neutrino effects. It constrains the number of thermalized neutrinos present at T ~ 0.1 — 1 MeV 
to be N u = 2.5 ± 0.7 [23]. Unfortunately, the uncertainty in N v is controversial so not strict 
bound on it can be derived. Here, we will simply take as a reference value the prediction of 
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Figure 8: The effective number of neutrino species as a function of Tr for different mixing angles. 
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Figure 9: il vs /^dm as a function of Tr for different mixing angles and m s = 1 keV. 

the standard cosmological model, N v = 3. Figure |8] shows N v as a function of Tr for different 
mixing angles. The variation with Tr is strong, going from N u ~ 3 — 4 for Tr > 7 MeV to 
iV„ ~ 0.3 for Tr = 1 MeV. The spread due to different mixing angles, on the other hand, is 
maximum (AN V ~ 1) at large Tr, and decreases for smaller Tr. Note that for sin 2 29 < 10~ 3 , 
N y is essentially insensitive to the presence of sterile neutrinos; it becomes a function only 
of Tr. As expected, the standard cosmological scenario is recovered at large Tr. In that 
region, if the mixing angle is large sin 2 29 ~ 0.1 all neutrinos -the three active plus the sterile- 
thermalize, yielding N u ~ 4. That is not necessarily the case for lower reheating temperatures, 
however. If Tr ~ 4 MeV, for instance, then N v ~ 3 for a sterile neutrino with sin 2 29 ~ 0.1; 
and the same N v can be obtained for sin 29 ~ 10~ 2 and T R = 5 MeV. Hence, LSND sterile 
neutrinos may still yield N v ~ 3, avoiding possible conflicts with big bang nucleosynthesis. 

The sterile neutrino relic density as a function of Tr is shown in figure |9| for different 
mixing angles and m s = 1 keV. Along the horizontal line, sterile neutrinos entirely account 
for the dark matter density of the Universe. The region above the horizontal line is therefore 
ruled out, whereas below it, v s only partially contribute to the dark matter density. Thus, 
in the region 3 MeV < Tr < 7 MeV and 10~ 3 > sin 2 29 > 10 -4 a sterile neutrino with 
m s = 1 keV may explain the dark matter. 

Because Q Ug scales linearly with m s , the results for a different value of m s can easily 
be obtained from the same figure. First notice from the figure that the sterile neutrino relic 
density also depends linearly on sin 2 29. So, another region where Q Vs = f^ m is m s = 10 keV, 
3 MeV <T R <7 MeV and 10" 4 > sin 2 29 > 1(T 5 . 

In the standard cosmological scenario, dark matter sterile neutrinos are produced at 
T ~ 150 MeV where collisions dominate the evolution of the neutrino system and matter 
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and thermal effects become relevant. As a result, the sterile neutrino relic density depends 
quadratically on m s and keV sterile neutrinos with sin 2 29 ~ 10~ 8 are required to account for 
the dark matter. In models with low reheating temperature, on the other hand, Q, Vs depends 
linearly on m s and much larger mixing angles are required to explain the dark matter. 

Cosmological and astrophysical observations can be used to constrain sterile neutrinos as 
dark matter candidates. The observed clustering on small scales of the Lyman-a forest, for 
instance, constrains the sterile neutrino mass from below. To obtain a limit on m s , the flux 
power spectrum of the Lyman-a forest must be carefully modeled using numerical simulations. 
The analysis presented in [|j and || respectively cite m s > 10 keV and m s > 14 keV as their 
limits, though a 30% discrepancy between them still exists. Such bounds, however, were 
obtained for sterile neutrinos produced in the standad cosmological model and do not direcly 
apply to the scenario we consider. That is why we will be mainly concerned with another 
bound, that derived from X-rays measurements. Sterile neutrinos may radiatively decay 
through v s — > v a + 7 producing a monoenergetic photon, J5~ = m s /2. X-ray observations 
may therefore be used to constrain or infer the mass of the sterile neutrino. In a recent analysis 
of the X-ray background from HEAO-1 and XMM-Newton, for example, the following limit 



relating sin 2 29, m s and fl l/s was found |11]. This bound is model independent, it applies to 
both the standard production mechanism and to the production in models with low reheating 
temperatures. 

In figure |l^ we display the sterile neutrino relic density as a function of sin 2 29 for different 
values of Tr and m s = 1 keV. The limit from X-rays, equation fl4.1|), is also shown and 
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Figure 11: The evolution of the sterile neutrino energy spectrum for Tr = 4 MeV, b = 10~ 3 and 
sin 2 26» = 1CT 8 . 

rules out the upper-right part of the figure. The different lines represent different reheating 
temperatures. Notice, for instance, that Tr = 4 MeV, £l Us = £ldm is n °t a viable point of 
the parameter space as it is incompatible with the X-rays limit. Indeed, sterile neutrinos can 
account for the dark matter only if Tr > 7 MeV. 

Turning this argument around we can also say that if dark matter consists of sterile 
neutrinos, they provide the strongest constraint on the reheating temperature. The present 
bound, in fact, gives Tr > 2 — 4 MeV and is based on the effect of active neutrinos on big 
bang nucleosynthesis. Dark matter sterile neutrinos might yield a more stringent constraint. 
Finally, notice that this bound on Tr was obtained for a sterile neutrino with m s = 1 keV 
but it only becomes stronger for larger masses. Dark matter sterile neutrinos, therefore, are 
useful probes of the early Universe. 

5. Sterile neutrino production in oscillations and decays 

The field (ft responsible for the reheating process may also have a direct decay mode into 
sterile neutrinos (4> — * v s v s ), opening an additional production mechanism for v s . As we will 
see, this mechanism significantly alters the predictions obtained in the previous section. In 
pEj], the production of sterile neutrinos in inflaton decays was investigated, but not in the 
context of low reheating temperatures. The main motivation to consider this mechanism 
is the conflict between the constraints from X-ray observations and those from small-scale 
structure that rule out the minimal production scenario for sterile neutrino dark matter. 
As mentioned in section 2, the decay eft — ► v s v s gives the following contribution to the 
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Figure 12: The sterile neutrino relic density 
as a function of sin 2 29 for T R = 4 MeV. The 
sterile neutrino mass is set to 1 keV and the 
curves correspond to two different values of b. 
The bound from X-rays observations is also 
shown. 



Figure 13: The sterile neutrino relic density 
as a function of sin 2 26 for T R = 4 MeV. The 
sterile neutrino mass is set to 10 keV and the 
curves correspond to two different values of b. 
The bound from X-rays observations is also 
shown. 



sterile neutrino collision integral 



2vr 2 

C^ Vs v B = b -^—j-^ T^dip - , (5.1) 



where b denotes the (ft branching ratio into sterile neutrinos, and m^, are respectively 
the (ft mass and number density. Being (ft non-relativistic, each v s is born with momentum 
p = m0/2, as enforced by the delta function. Due to this new contribution, f Vg will now 
depend not only on Tr, m s , and sin 2 26 but also on b and m^. To keep things simple we will 
set rrifj, = 100 MeV and study the dependence of the different observables with b. 



Figure 11 displays the evolution of the sterile neutrino energy spectrum for Tr = 4 MeV, 
b = 10~ 3 , and sin 2 26 = 10 8 . Each line corresponds to a different temperature. It is not 
difficult to decipher what is going on. Whenever a (ft decays, a peak at p = m^/2 in f Va 
is generated. But not all eft's decay at the same time. And the momentum of the sterile 
neutrinos produced in earlier decays is redshifted when later decays occur. That is why, at 
any given temperature, the resulting spectrum has a drastic jump at p ~ m ( f ) /2, with all 
the neutrinos produced before (in decays) lying at smaller momenta. As we approach the 
radiation dominated epoch, the redshift essentially ceases and only residual decays modify 
the spectrum at large p/T 7 . At the end, no traces of the discontinuity at p = m^/2 are left 
in the primordial spectrum. 

The sterile neutrino relic density is shown in figure 12 as a function of sin 2 26. For that 



figure Tr = 4 MeV, m s = 1 keV and the two curves correspond to b = 10~ 2 and b = 10~ 3 . 
The solid line is the X-ray constraint obtained from equation ( |4.l[) . The relic density behaves 
in a similar way for the different values of b. At large mixing angles, the production of sterile 
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neutrinos is dominated by oscillations and independent of b. That is the case we dealt with 
in the previous section. At smaller mixing angles, we encounter an intermediate region where 
both production mechanisms are relevant and the relic density depends on b and sin 2 29. 
Finally, at even smaller mixing angles, sterile neutrinos are produced dominantly in <fi decays 
and therefore the relic density does not depend on sin 2 29, as signaled by the horizontal lines 
observed in figure [l^. In that region the sterile neutrino relic density is simply proportional to 
b. If sterile neutrinos account for the dark matter, Vt Va = fi<j m , the X-rays constraint requires 
a small mixing angle, sin 2 29 < 10~ 4 . 

New viable regions, where the sterile neutrino is produced in <fi decays and makes up the 



dark matter of the Universe, can be read off figures 12 and |13j. For instance, a m s = 1 keV 
sterile neutrino with sin 2 29 < 10 -4 will be a good dark matter candidate for Tr ~ 4 MeV 
and 10 -3 < b < 10 -2 . For decay-dominated production, Q Us is simply proportional to Tr, 

VL Vs oc bm s T R . (5.2) 



Using this equation in conjuntion with figures |12| and |13, additional allowed regions can be 
found. 



Figure 13 is analogous to figure 12 but for a larger value of the sterile neutrino mass, 
m s = 10 keV. The two curves correspond to b = 10 -3 and b = 10~ 4 . Owing to the increase 
in m s , the X-ray limit becomes much stronger than in figure [l^. Indeed, it constrains dark 
matter sterile neutrinos to have a very small mixing angle, sin 2 29 < 10 -9 . 

In the standard production mechanism, such small mixing angles are not allowed as they 
yield a too small sterile neutrino relic density, Q Us oc sin 2 29. For sterile neutrinos originating 
in (p decays, on the contrary, the production mechanism and the radiative decay are controlled 
by two different parameters. In fact, Q Ua oc b whereas T(u s — ► v a + 7) oc sin 2 29. Thus, no 
matter how small sin 2 29 -and consequently T{v s — > v a + 7)- is, it is still possible to find 
appropriate values of b, Tr and m s such that £l Ug = £ldm- In other words, for b ^ the 
X-rays limit can always be satisfied. 



6. Conclusions 

We numerically studied the production of sterile neutrinos in models with low reheating 
temperatures. Two production mechanisms for the sterile neutrinos were taken into account: 
active-sterile neutrino oscillations {y a «-> u s ) and the direct decay of the scalar field (0 — > v s v s )- 
Several observables, including f Ug , N u , and Q Va , were computed for different sets of reheating 
and mixing parameters. We showed that in these models, LSND sterile neutrinos may still give 
N v ~ 3 -avoiding problems with big bang nucleosynthesis- and that keV sterile neutrinos 
may account for the dark matter of the Universe. Dark matter sterile neutrinos produced 
in oscillations were found to be effective probes of the early Universe, as they constrain the 
reheating temperature to be rather large, Tr > 7 MeV. Finally, we showed that sterile 
neutrinos originating in decays may explain the dark matter and satisfy the bounds from 
X-ray observations. 
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